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We consider a crowd of N persons trying to exit some area trough a small exit. The probability is calculated that an individual 
is able to withdraw from the crowd under one's own steam. The problem is simulated within the generalized force model (D. 
Helbing et al., Nature 407 (2000) 487), and all model parameters are taken from this paper. The results indicate, that in a crowd 
of 150 persons, this probability is not greater than ten percent. We also evaluate the number of helpers necessary to get the above 
O-i probability of fifty percent. 
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1. Introduction 



When applications of natural sciences to human beings 
are considered, the problem of unpredictability of human 
mind is an eternal motif. In psychology, neurophysiological 
arguments overlap with philosophical ones [1, 2]. In sociology, 
the empiricistic foundations of natural sciences have been 
^ jindermined by interpretative sociology of Weber and Simmel 
^tOll; 'Determinism is dead' is a catchphrase of today 0], 
Against this background, the modeling of crowd dynamics - 
, which is our aim here - could seem to be an extreme example 
of a mechanicistic reductionism. On the other hand, prediction 
y—i is an ultimate aim of all sciences; here the famous statement 
on scientific aim 'Savoir pour prevoir et prevoir pour pouvoir' 
by Auguste Comte remains valid and desirable [5J]. In fact, 
a crowd can be seen as a many-body system with local 
interactions; in such systems, statistical laws should allow for 
some predictions. Once a comparison with experimental data 
became possible the critique from the hermeneutically 

oriented audience is less convincing. 
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Modeling of the crowd dynamics is known for more than 
50 years ]9[]. The methods have been advanced much in 
90's by Dirk Helbing and cooperators; reviews and lists of 
references can be found in ifiol [Til f*7tl _ Among the methods, 
the social force model seems to be most realistic [12]. This 
is a set of differential equations of motion, where positions 
and velocities of pedestrians are time-dependent variables; 
the approach is equivalent to the molecular dynamics, where 
human desires are encoded in the form of social forces. Within 
this model, a number of problems have been addressed, as 
lane formation, strip formation, turbulent waves, herding and 
bottlenecks Hill . In this paper we address more directly to 
the problem of the phase of a simulated crowd. Namely, we 
ask for the conditions when an individual or a small group can 
change their position with respect to their neighbours in the 
crowd. If they are stuck, we refer to a clogged phase. This 



criterion is a direct analogy to the Monte Carlo simulations of 
the crystallization phase transition of hard spheres lfl3l [l4Tl . 
We note that although the dynamics of hard spheres is much 
simpler than the crowd dynamics, systematic numerical stud- 
ies of the finite size effect have been possible only recently 11511 . 



The scenario to be simulated here within the social force 
model is as follows. Pedestrians numbered by i = l,...,N 
are going to leave a room through a small exit. During this 
process, we monitor the sum S , of mechanical compressive 
forces acting on each individual. Once this sum exceeds some 
prescribed value S c for any individual j, the direction of the 
desired motion of this individual is reverted from the vector 
towards the exit to the opposite. Further, a number K of 
pedestrians who are nearest neighbours of y'-th one decide to 
accompany her/him. Then their directions of desired motion 
are set equal to the direction of j-th individual, and their social 
forces towards j-th individual, initially repulsive, change signs. 
The outcome of the simulation is the probability P, that the 
crowd throws j out through the exit, despite her/his struggling 
to withdraw from the crowd. As K pedestrians help j, P is 
expected to decrease with K. 

Among the problems considered by other authors, this 
scenario is somewhat similar to the bi-directional flows in 
the bottleneck problem ifTHl . where two flows of pedestrians 
walking in opposite directions met at a narrowing of the path. 
Also, the effect of clogging has been demonstrated in Illa. | l7ll 
as long time periods when nobody crosses an exit. In 11811 . 
a picture is drawn (Figs. 4 and 5) where a narrow queue 
of pedestrians is formed through the crowd of individuals. 
The effect of the crowd solidification as opposed to the lane 
formation was discussed in |T^[. In this paper, the phase of 
coherent motion in lanes was destroyed by an added noise. 
In i20ll . the authors defined clusters of individuals as sets of 
those who interacted via physical forces. Above some optimal 
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value of the desired velocity, the distribution of size of these 
clusters was found to be strongly modified. In our earlier paper 
1 2JJ] we discussed a similar scenario, where pedestrians around 
make more place to a handicapped person, for she/he can 
leave the room. In the text presented here, a single individual 
or a small group intends to move with respect to the crowd; 
then, the crowd itself is considered as a medium for individual 
pedestrians or their small groups. Up to our knowledge, this 
issue has not been discussed yet. 

The paper is organized as follows. In the next section the 
social force model is explained in more detail. We adopt the 
formulation and the values of the parameters used in [ 12]. Sec- 
tion 3 provides the numerical results on the probability P of 
being stuck as dependent on the model parameters N, K and 
S c . Last section is devoted to concluding remarks. 



2. The model 

The model equations of motion are adopted from the gener- 
alized force model Il8lll9lll2ll . Such an equation for a person 
of mass m is as follows 



(hi v(r,) - Vj 

m — = m 

dt t 



(i) 



where the first term on right hand side is the tentative accel- 
eration of a person i who intends to have the velocity v(r,), 
dependent on the coordinates r,-; as a rule, the vector v points 
to the exit center (large distance from the person to the exit) 
or to the closest point of the exit (small distance). In our 
simulations, the absolute value of the desired velocity v(r,) 
is 3 m/s; this is purposefully higher than the optimal value 
1. 375 m/s l20ll . Further, r is the characteristic time of this 
acceleration, v,- is the actual velocity of i-th person, fy is the 
force exerted on i-th person by y'-th person, and f,-u/ is the force 
exerted on i-th person by a wall W. The force fy contains three 
components; 'psychological' interaction which describes the 
tendency of i and j to keep distance between each other, and 
two physical interactions between their bodies: radial force 
and slide friction. The psychological interaction is equal to 
A,exp((2/? - ||r,- - r/||)/B), where r, is the position of i-th 
person, R is the mean 'radius' of the vertical projection of the 
human body. This psychological part of fy will be modified 
in our simulation. The radial physical force is equal to 



kg(2R- 



where g(x) = x if x > 0, g(x) - elsewhere, 



and A: is a constant. The physical friction is assumed to be 
Kg(2R - ||r,- - rj )((Vj - V,) ■ tjj)tij, where t,j is the unit vector of 
tangential direction to the body surfaces. The same expressions 
of the physical forces are used to describe the body-wall 
interaction. The instant values of the velocities V; allow to 
update the positions r,- as well. The parameters of the system 
are adopted from HI. Namely, A = 2000 N, B - 0.08 m, t = 
0.5 s, k =1.2 xlO 5 kg/s 2 , K/k = 2, R = 0.3 m and m = 75 kg. 
The room is 30 x 20m, with the exit of width of Im in the mid- 
dle of the shorter wall. The initial positions of the agents are 



Figure 1 : A spatial configuration of individuals near the exit. X is marked with 
square (blue online), the six helpers are marked with triangles (blue online), the 
others are with circles (red online). In this example, the help is successful. 
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Figure 2: The probability P of being thrown through the exit as dependent on 
the number of persons N in the crowd, for different values of the threshold S c . 
These data are obtained for K = (no helpers). 



selected randomly, but the initial overlaps of them are excluded. 

To determine if a person is stuck or not, a numerical exper- 
iment is performed as follows. The sum 5/ of compressive 
mechanical forces acting on each individual i is registered 
during the motion. Once for some individual X the sum S x 
exceeds some threshold value S c prescribed at the beginning 
of the simulation, the desired direction of motion for this 
individual is reverted. Now this direction is not to the exit 
but the opposite one. We assume that all individuals act with 
the same force. We are also interested if a collective action 
of neighboring individuals could change the outcome of the 
experiment. To check this, we change also the desired direction 
of K individuals, who are closest to X when the threshold value 
is exceeded. Their desired direction of motion is now equal to 
the desired direction of X. Also, the repulsive psychological 
forces between X and his neighbours change sign to be 
attractive. Now the group of K + 1 individuals tries to evade the 
exit, as if they tried to help a victim of the interpersonal forces 
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in the crowd. In Fig. 1, an example is shown for K — 6, where 
the help is successful. In both experiments, if X crosses the 
exit despite this change of her/his intention, we call the crowd 
'jammed'. 



3. Results 

The outcome is the probability P that X is thrown out 
through the exit against her/his will. In Fig. 2 this probability 
is shown for K = (no helpers) against the crowd size N, for 
various values of the threshold value S c . Each point on these 
results is an average over 70 samples. As we see, the results 
only weakly depend on S c . On the contrary, the crowd size N 
is decisive. As we see, the probability P increases with N from 
less than 0.4 for N = 50 to about 1 .0 for N = 250. 

The simulations are repeated in the presence of K helpers, 
for K between 1 and 10, S c — 150 N. The results are shown 
in Fig. 3. Each point is an average over 200 samples, and the 
error bars are the differences between the averages over first 
and second hundred of samples. Here again, N is relevant, 
but the crowd size can be to some extent neutralized by the 
number of helpers. For example, P close to 0.5 can be achieved 
in a crowd of N = 100 persons with K about 2 helpers, in a 
crowd of N - 200 persons with K about 5 helpers and so on. In 
simulations, we observed many times that the group of helpers 
become dispersed, as in Fig. 1. In the presented results this 
dispersion is neglected; this means that the effective number of 
helpers is perhaps smaller. 
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Figure 4: The probability P of being thrown through the exit against the number 
of helpers K, for different values of the threshold force 5 c . Here, the crowd size 
is N = 150. 



were N = 150, K = and S c = 150 N. The obtained probability 
P, which is 1.0 for r = 0.5, is not less than 0.95 for r as large 
as 5.0. We deduce that the variation of r is not relevant. 



4. Discussion 

Our numerical results indicate that once the crowd size N 
exceeds 150-200 persons, it is unlikely that a single individual 
can withdraw under one's own steam. Then, any large gath- 
ering of people should be treated as a potentially dangerous 
medium. This conclusion is supported by the historical data on 
crowd disasters [22]. 
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Figure 3: The probability P of being thrown through the exit against the number 
of helpers K, for different values of the crowd size N. Here, the threshold force 
S c = 150 N. 

The simulations with helpers are repeated also for N = 150 
and different S c . These results are shown in Fig. 4. Each point 
is an average over 100 samples. These results confirm, that the 
threshold value S c does not influence much the probability P. 
We have performed also a similar experiment with the variation 
of the time parameter r. Other parameters of this experiment 



In the clogged phase, the only chance to leave the crowd is 
to mobilize a group of helpers nearby. We note that once this 
group is dispersed, each separate person is helpless in the same 
way. What does matter in these conditions is the communica- 
tion between people. The question arises, how many helpers 
must be found to have a chance of 50 percent to withdraw from 
the crowd. Let us denote this number as K^o(N). Because of 
the complexity of the problem, the accuracy of our results al- 
lows to classify them as semi-quantitative only. We made an 
attempt to fit log(Kso) against log(N) to obtain the exponent /3 
in the tentative scaling relation K$o °c N 13 . The result is that 
P =1.88 ± 0.05. We feel entitled to claim that /3 is larger than 
1.0. However, as cannot be greater than N, this behaviour 
must end with some crossover for larger N. A more quanti- 
tative evaluation of the function shape of K^oiN) needs much 
more computational power. 

Summarizing, in a crowd of some hundreds of people the 
difference between an autonomous human being and a piece 
of passive body is less than we would like to admit. Although 
at the exit the individuals leave the room one by one, in the 
middle of the crowd their mutual positions cannot be changed. 
In these conditions, the unpredictability of the human mind 
does not influence the trajectory of its owner. The obtained 



3 



20 




1 1 1 ' ' 1 

100 150 200 250 300 

N 

Figure 5: The probability P of being thrown through the exit against the number 
of helpers K, for different values of the threshold force S c . Here, the crowd size 
isN= 150. 

results should be helpful to evaluate human resources which 
are needed to tackle emergency situations in large gatherings 
of people. 
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